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ABSTRACT 


A  dynamical  analysis  of  the  M60AIE2  Tank  main  gun's  elevation  stabiliza¬ 
tion  system  is  performed.  The  equations  describing  the  gun,  hull,  and  turret 
dynamics  are  derived  and  computer  diagraming  is  shown  for  simulation  purposes. 
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1.  Introduction 

The  M60AIE2  Tank  stabilization  gun  control  system  was  designed  to 
stabilize  inertially  the  turret  and  cupola  in  traverse  and  the  main  gun  and 
commander's  machine  gun  in  elevation.  These  are  depicted  in  Figure  1.  The 
four  stabilization  loops  including  the  turret  traverse,  cupola  traverse,  main 
gun  elevation,  and  cupola  gun  elevation  are  similar.  Each  uses  electro  hydrau¬ 
lic  actuators  to  drive  the  load  based  on  information  derived  from  the  gunner's 
controls,  commander's  controls,  or  rate  gyro  sensors  in  the  cupola,  hull,  and 
turret.  Because  of  a  variety  of  contributing  factors,  the  stabilization  system, 
as  originally  designed  and  installed  by  the  contractors,  had  performance 
deficiencies. 

The  task  of  the  Army  Inertial  Guidance  and  Control  Laboratory  and 
Center  was  to  investigate  the  system  dynamics  and  to  develop  fixed  gain  compen¬ 
sation  networks  for  the  stabilization  loops.  The  initial  efforts  were  directed  at 
investigations  in  the  azimuth  plane  only.  In  it,  the  turret  and  hull  assemblies 
were  analyzed,  and  the  model  was  synthesized  on  two  EAI-22IR  analog  computers 
for  further  study. 

At  the  request  of  the  Program  Manager,  a  scope  of  work  for  a  follow-on 
effort  was  submitted.  In  essence,  the  task  was  to  cover  the  development  of 
improved  fixed  gain  compensation  networks  for  the  main  gun's  elevation  loop. 


2.  Purpose 

It  is  the  purpose  of  this  report  to  show  the  development  of  the  dynam¬ 
ical  equations  of  motion  of  the  main  gun,  turret,  and  hull  in  the  elevation  plane. 

In  addition,  the  synthesis  of  the  mathematical  model  will  be  displayed  for  the 
analog  computer  simulation. 


3.  Main  Gun  Dynamics  in  the  Elevation  Plane 

The  physical  situation  of  the  main  gun  within  the  tank  turret  is 
shown  in  Figure  2.  The  geometry  and  free  body  diagram  of  the  forces  acting 
on  the  gun  are  shown  in  Figure  3. 

The  dynamics  are  described  by  summing  all  the  external  torques  on  the 
main  gun  and  setting  them  equal  to  the  rate  of  change  of  angular  momentum  of 
the  gun: 
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FIGURE  1.  THE  M60AIE2  TANK  STABILIZATION  GUN  CONTROL  SYSTEM 
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First  derivative  of  momentum  with  respect  to  the 
mertial  haais  )  . 
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From  the  right  side  of 


External  torques  acting  on  the  gun  summed  about 
its  mass  center. 

equation  ( 1 ) ,  the  momentum  is  given  by: 


HG  =  IG  •  -G'!  (2) 

where 

1^  ^  Inertia  tensor  of  the  gun 

,  Q  j 

w  ^  Angular  velocity  of  the  gun  with  respect  to  the  inertial  basis. 

Q— I 

Now  oT  has  two  components  in  this  analysis;  i.e., 


where 


—  G-I  — G-H  —H-I 

w  =  cu  +  cu 


(3) 


^  Q^H  A 

oT  =  Angular  velocity  of  the  gun  with  respect  to  the  hull 

*  =  Angular  velocity  of  the  hull  with  respect  to  the  inertial  basis. 

Since  the  hull  is  also  described  within  a  rotating  reference  frame,  there  is  a 
Coriolis  force  to  be  considered. 
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(4) 


—  G-I 

oj  x 


where 


S-  d  h" 

HG  dt 


Fi'  .it  derivative  of  momentum  with  respect  to  the  gun 
i*  ->is  (XG  ,  Yg) 


Substituting  equations  ,  2)  and  { 3)  into  { 4)  and  performing  the  indicated 
differentiation  giver  '  ie  following: 


I 

h\ 


d_ 

dt 


?  —  G-H  — H-I 

I  •  U)  +  cu 

(j 


/  —  G-H  ^  .  I-  G-H  _H-l\ 

+  J  x  IG  +  u  j 

(5) 


G 

=  Z. 


( —G-H  —H-I  \ 

•  la  +  u>  j 

r  f-G-H  S-K-d 
lG  *  \w  +  u  / 

/—G-H  —H-I  \  ^ 

x  )*  *G  ’ 


—G-H  —G-I 

c o  *■  cu 


'  (6) 


The  first  term  on  the  right  of  the  equality  sign  in  equation  ( 6)  goes  to  zero 
G 

because  the  term  IG  ,  the  derivative  of  the  inertia  of  the  gun,  is  zero  in  the 
gun  bar  is. 
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Now  from  equation  (1) ,  develop  ihe  external  force  acting  on  the  gun. 


7.  Mi 

—  cm 
i=l  f 


a 

_  V 

LJ 

i=l 


ML 


gun 


(a+0) 


(8) 


where 


n 

Yj  Mi^,  =  Moments  of  the  gun  about  the  trunnion 

i=l  gun 


MQif  (a  +0)  & 


The  inertia  force  of  the  mass  center  displaced  by 
distance  i  i  from  the  trunnion. 


n 


=  MGgf3  + 


F  cos  0  Li 

Pi 


(9) 


But 


13  =  cos  a 
=  (i2  - ij)  cos  or  . 

Thus 

n 

y  Mi  =  M  gi  j  cos  a  +  F  cos  /3  (t.  2  - !  cos  a 
.  .  t  Cj  A 

1=1  gun 

+  MGif(a+9)  .  (10) 

Combining  equations  (7)  and  (10)  gives,  on  the  right,  the  invariant  vector  form 
of  the  gun  dynamics. 
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Equation  ( ll)  is  still  not  usable  for  computer  synthesis  because  the  right  hand 
side  is  still  in  vector  form.  It  would  be  desirable  to  have  the  vector  terms 
reduced  to  scalar  format  in  the  gun  basis.  This  is  done  next.  The  subscript 
on  the  inertia  terra  will  be  dropped  to  avoid  confusion  with  the  notation  of  a 
matrix  in  the  gun  basis. 


But  to  be  compatible  with  the  left  side  of  equation  (li) ,  the  vector-matrix 
form  of  equation  (13)  must  be  transformed  from  the  H-basis  to  the  G-basis: 


=  C 


Thus, 


cos  y  sin  y  0 

C  =  -sin  y  cos  y  0 

%  n  n  i 


cos  y  sin  y  0  0 

-sin  y  cos  y  0  0 


This  result  is  not  as  striking  as  it  first  may  appear.  It  says,  in  effect,  that 
there  is  no  difference  between  the  H-basis  and  the  G-basis  for  the  angular 
rate  shown.  This  of  course  is  not  true  since  the  two  bases  are  not  equivalent. 
What  is  actually  seen  is  that  the  concept  of  a  "pseudo-vector, "  namely  the 
angular  rate,  is  in  a  third  dimension,  but  the  oniy  concern  here  is  in  the  X,  Y 
plane  for  this  analysis.  Thus,  that  particular  pseudo-vector  is  not  transformed 
only  for  this  special  case. 


It  remains  to  reduce  the  Coriolis  term  to  its  scalar  equivalent  in  the 
G-hasis: 


The  right  side  is  expanded  as  follows: 


=  ~tl  (a  +6)  +  jl  (a+0)  ,  (19) 
yy  aa 

and  the  vector  inner  product  is  given  by 


0 

0=0  .  (20) 
a  +  0  _ 

Once  again  the  problem  of  a  two-dimensional  simulation  while  working 
with  a  three-dimensional  real  model  is  evidenced.  Equation  (20)  says  that 
there  is  no  Coriolis  term  contribution.  The  purpose  here  was  to  verify  a 
relation  that  is  well  known.  Cannon  [1]  says, 

For  plane  (two-dimensional)  motion  a  major  simplifica¬ 
tion  occurs  in  the  equations  of  motion  because  we  are  concerned 
only  with  axes  of  rotation  and  angular-momentum  vectors  in  a 
single  direction,  namely,  perpendicular  to  the  plane  of  motion; 

.L 

in  this  case,  of  the  36  possible  terms  in  H  ,  only  one  remains! 
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To  wit: 


I 

H  =  T  J  h  (21) 

2  7.  Z 


1^ 

in  which  H  is  the  rate  of  change  of  angular  momentum,  which 

is  now  perpendicular  to  the  plane  of  motion. 

An  even  more  graphical  representation  of  the  results  of  equation  ( 20)  is 
obtained  if  one  uses  the  notation  of  Cannon  on  equation  ( 17) , 

u  x  I  •  iu  =  uf  x  J1  •  u>t  =  0  .  (22) 

2  Z  Z 

The  strong  evidence  of  the  triple  scalar  product  of  collinear  vectors  being 
zero  is  proof  enough. 

Finally,  the  scalar  form  of  equation  (11) ,  in  the  G-basis,  is  given  by 

M  gfj  cos  a  +  F.  cos  /3  (i2  —  i j)  cos  a  +  Ml](a  +  0)  =  J_ ( a  +  0)  . 

G  G  G  (23) 

Another  dilemma,  which  is  not  apparent  from  equation  (23) ,  is  the 
source  of  the  information  to  obtain  a  and  9  .  As  previously  mentioned,  ex  is 
the  angle  between  the  gun's  x  reference  axis  and  the  hull's  x  reference  axis, 
and  0  is  the  angle  between  the  hull's  x  reference  axis  and  the  inertial  x  reference 
axis.  Obtaining  0  is  straightforward.  It  is  merely  picked  off  of  the  turret 
elevation  hull  gyro.  But,  a  is  nor  measurable  directly,  since  the  main  gun 
reference  gyro  also  measures  rate  with  respect  to  inertial  space.  However, 
this  can  be  alleviated  with  the  following  consideration: 


G-I  H-I  G-I  I-H  G-H 

a>  -  co  =  co  +  co  -co  =  a 

G  G  G  G  G 


(24) 


Thus,  the  angular  rate  of  the  gun  with  respect  to  the  hull  can  be  instrumented 
by  taking  the  difference  in  the  outputs  of  the  two  reference  gyros. 


4.  Vehicle  Dynamics  in  the  Elevation  Plane 

The  turret,  cupola,  and  hull  are  all  considered  to  be  a  single  rigid 
body  in  this  analysis.  This  is  not  true  actually,  but  the  cupola's  mass  and 
moment  of  inertia  are  so  small  compared  to  the  overall  tank  that  it  can  be 
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included  without  consideration  o£  its  change  for  various  hull  elevation  angles,  9  . 
There  is  also  a  negligible  effect  between  the  turret  and  the  hull  known  as 
"dishpanning. "  This  is  of  little  consequence  with  respect  to  the  degree  of 
accuracy  of  this  analysis.  The  geometry  of  the  tank  vehicle  is  shown  in  Figure 
4.  The  tracks  and  suspension  are  considered  as  a  pair  of  spring  and  damper 
combinations.  This  allows  two  degrees  of  freedom  in  this  plane  if  the  forward 
velocity  is  not  considered.  Figure  5  shows  the  free  body  diagram  of  the  forces 
on  the  tank. 


FIGURE  4.  GEOMETRY  OF  THE  COMPLETE  TANK  WITH  SUSPENSION 
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FIGURE  5.  FREE  BODY  FORCE  DIAGRAM  OF  THE  TANK  AND  SUSPENSION 

By  the  principle  of  D'Alembert,  the  summation  of  the  forces  in  the 
Y  direction  are  given  by: 


£f  ,  ~  F,  -  F,  -  F,  -  F,  -  f.  =  0 


From  geometry  in  Figure  4, 


Ft  =  KjY, 


F2  =  BjYj 


F3  =  K2Y2 


f4  =  B2Y2 


fi  =  V  • 


Yj  =  Y  -  fj'  sin  0 


Y2  =  Y  +  (£2  -  f !  )  sin  0 


Thus, 
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Finally,  simplifying  and  rearranging  yields, 

JJ  +  3cos20(B1i1,2-B2i-22  -B2/1!2+2B2f1'f2') 

+  cos  9  sin  0  (Kj/j'2  -  K2i2’2  -  K2f,’2  +  2K2fi,!2’) 

+  Y  cos  -  B2i2  +  B2ij')  +  Y  cose  (K *!/  -  K2i2  -  Kjf,’)  =  0  . 

(40) 


This  system  is  said  to  be  elastically  coupled  because  of  the  unsymmetric 
location  of  the  mass  center  with  respect  to  the  suspension  system. 

It  should  be  noted  that  equations  ( 34)  and  ( 40)  are  exactly  accurate  only 
for  the  case  in  which  the  tank  has  no  roll  with  respect  to  iuertail  basis  asso¬ 
ciated  with  its  motion.  This  is  reasonable  only  if  it  is  assumed  to  be  traveling 
on  a  model  of  level  terrain  with  only  pitch  undulations. 


5.  Actuator  and  Main  Gun  Dynamics 

For  the  case  of  no  compliance  in  the  rod  or  in  the  fluid,  Figure  6 

shows  that 


Y 


l 


=  Y2  =  Y  . 


(41) 
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Thus,  fiuid  flow  into  the  valve  gives  a  force  equation 

Q  =  Y,A  ,  (42) 

where 

Q  =  fluid  flow  rate 

Yj  =  rate  of  displacement  of  the  valve 

A  =  area  of  the  valve  face. 

Similarly,  with  no  compliance,  and  a  small  angle  approximation, 

Y  2  =  la  .  (43) 

Combining  equations  (43)  and  (44)  yields 

Q  =  f  a  A  (44) 

The  general  hydraulic  equation  including  losses  given  by  Wroble  [  2]  and 
Cannon  [1]  is  modified  to  give 

Q  =  KY0^P  -P.  -AP  .  (45) 

S  JLi 

Thus, 

Iff  A  =  KY0n/p  -  P  -  AP  .  (46) 

S  JL 

and 

a  =  K*  Yq'J P  -  P  -  AP  ,  (47) 

S 

where 


K*  -  JL 

K  i  A 

P  =  supply  pressure  to  actuator 
s 

P  =  losses  in  pressure 
Li 

A  P  =  change  in  pressure  measured  across  the  valve  . 
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6.  Mode!  of  the  Actuator  -  Equilibrotor 


The  equilibrator  and  actuator  combination  is  shown  in  Figure  7. 


FIGURE  7.  EQUILIBRATOR  AND  ACTUATOR  COMBINATION  I 

i 

The  tree  body  force  diagrams  are  shown  in  Figure  8. 


FIGURE  8.  FREE  BODY  FORCE  DIAGRAM  OF  ACTUATOR 
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The  problem  is  to  find  the  force  F  =  F  as  a  function  of  known  quantities 

1 1  K1 

including  F,  Y.  and  a.  As  shown  in  Figure  8,  parts  (b)  and  (c),  the  springs 
and  piston  each  have  no  appreciable  mass;  i.e. ,  m  =  0  . 

From  the  free  body  diagrams, 

F  =  F  (48) 

Yi  K1 


Thus 


F  =  F  +  Fc 
Ki  K1 


F  t  Fc  =  F 
K2 


F  =  F  +  Fc  =  F 
Y1  K2 


FK1  =  K’Y^ 

FK2  =  K^Y-Y,) 

Fc  =  Fc  Sgn(Y  -  Yt) 
max  °  1 


Now,  Y  is  an  obtainable  quantity,  while  Yj  is  not.  The  expression  for  Y  is 
obtained  from  solution  of  equations  (34)  and  (41).  Therefore,  Yt  must  be 
replaced  by  a  known  value;  from  equation  ( 53) , 


Y,  = 


Thus, 


kKv-fkj) 


Fc  =  Fo  Sgn  Y  -  ~ 

max  &  dt  \  Kj 
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Substituting  equations  (57)  and  (58)  into  (51)  yields: 


4  -  V)  ■ 


Fc  Sgn  (y  -  j-  --  F 

max  \  dt  Kj  J 


Solving  for  F  : 

K.1 


.  „  •  FK1  „  _  *  d  { FK1  \ 

K2 Y  -  K2“kT  +  i?CmaxSsnLY'drV'ir/J  =F 


dt  \  k Y) 


KjL  ,  _  _  •  d  {  FK1 

-  -  -rrf-  {  F  -  K2  Y  -  Fc  Sgn  Y  -  —  l  —rr~ 

!  *•  mav  °  dt  \  K* 


K1  K2 


Now  the  expression,  F,  is  given  by  , 


F  =  1.885^-P^  -  0.  587(1740  +  7.5a)  +  20(0.783)  1  (61) 


Thus, 


=  -  -§f-  (  1 . 885^P_  -  -  0.587(1740  +  7.5a)  +  20(0  783) 

K1  K2  \  U  b/ 

/F  \  | 

-  K2  Y  -  Fc  Sgn  Y  -  {~^}  }  .  ( 

1  max  \  Kj  /  j 


Differentiating  equation  (63)  with  respect  to  time  gives  the  expression  for 
inclusion  in  the  signum  function  of  the  coulomb  friction  term: 


d  ("sf)  =  '  kIkT  ['°-587(7-5“)  "k2^] 


dt  \  KJ 


so  that  finally, 


F  =  F 
A  K1 


--§f-  1 1 . 885/jP  -  P_\  -  0.587(1740  +  7.5a)  +  20(0.783) 

K2  1  \  >' 

K2 Y  "  Fcmax  Sgn  [Y  “  (°-  587(  7l  5c*>  +  K2 


’Equation  (61)  was  derived  by  T.  G.  Wetheral. 
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Now,  equation  (65)  is  exactly  the  expression  needed  to  complete  the  moment 

equation  of  the  gun  given  by  equation  ( 23) ,  where  the  actuator  force  F  in 

that  derivation  is  the  same  force  as  F  here.  ^ 

K1 

The  alternative  model  of  the  Actuator-Equilibrator  is  shown  in  Figure  9. 


FIGURE  9.  ALTERNATIVE  ACTUATOR  AND  EQUILIBRATOR  COMBINATION 
The  free  body  force  diagram  is  shown  in  Figure  10. 


FIGURE  10.  FREE  BODY  FORCE  DIAGRAM  OF  ACTUATOR 


From  the  hydraulic  equation  this  relationship  is  given 

Yi  -  Y0  =  KX\/f  -  F  -  B  Sgn(Y,  -  Y0>  (65) 

3  c 

and, 

F  =  ( Yj  -  Y2)  Kf  +  ( Yj  -  Yj)  Bf  .  (66) 

This  can  be  visualized  in  Figure  11. 2 


FIGURE  11.  ANALOG  DIAGRAM  OF  ALTERNATIVE  SCHEME 


Now 


Q  =  11  gal/min  x  x 


231  in.3min 


=  42.35 


in. 

sec 


SO  gal  sec 

at  the  third  stage  displacement  of  the  servo  spool  of  X0  =  0. 125  inch.  Since 


(67) 


Q  in.3/sec  42.35  in.  J  .  . 

Y,  =  .  — T— T—  =  tt7 - =  1.872  ft/sec  . 


A  in. 


1.885  sec 


(68) 


^This  derivation  was  done  by  W.  C.  Jordan. 
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The  supply  pressure  is  given. 


F  =  2000  lb/in.z  x  ir885  in.2  =  3770  1b 


fl  -  0 


since  Q  is  the  no  load  flow.  Also,  the  hydraulic  equation  is 

Q  =  KXon/Y  -  F  =  KXosfT  ; 

”  S  L  —  s 


Yj  =  K’x  0.125  in.  x  s!  3770  lb 


K'  =  0.2439  - ^ - 

in.  -sec 41b 


7.  Anoiog  Computer  Synthesis  of  Moin  Gun  Dynamics 

The  dynamics  of  the  main  gun,  given  by  equation  (23) ,  are  repeated 
here  as  follows; 

••  •»  1  «■ 

(a  +0)  =  (a  +9)  MGf?+ cos  +  FA  cos^  (i2 -i,)j  (72) 

This  is  represented  by  the  circuit  in  Figure  12. 


8.  Analog  Computer  Synthesis  of  Tonic  Dynamics  in  the  Pitch  Plane 

The  linear  motion  of  the  tank  in  the  pitch  plane  (with  coupling)  is 
given  by  equation  (34) ,  rewritten  here 


-y  -  Tr-[uB, 

H 


+  B2)  Y  +  (Kj  +  K2)  Y  -  (Kji/  +  K2X2’  -  K2i/)sin  0 


(Bjii'  +  B2ii  -  B2i2')  cos  99]  . 
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(73) 


Jc 


i?) 


G  1 


This  equation  is  programmed  as  shown  in  Figure  13. 

The  angular  motion  of  the  tank  is  also  coupled  into  the  linear  motion  in 
the  pitch  plane.  The  expression,  rewritten  from  equation  (41),  is 


-0  =  ~~~  f[  0  cos  0  (  Bj/j'2  -  B2f  2  2  ~  B2i  j'2  +  2B2f  i  f  2 ) 

JH  L 

+  cos  0  sin  0(K1l/2  -  K 2f22  -  K +  2K2fii2') 

+  Y  COS  *"  ^>2^-2  ®2^l) 

+  Y  cos  0(K2fj  -  K2f2  _  Ktij')]j=  0  .  (74) 

This  equation  is  programmed  as  shown  in  Figure  14. 

These  analog  diagrams  representing  the  gun  and  hull  dynamics  will  be 
incorporated  into  the  system  simulation  shown  in  Figure  15. 
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<*/,2*  K2i'22.  K2i'22  +  2K2i'  1'2) 

FIGURE  14.  ANALOG  DIAGRAM  OF  TANK'S  COUPLED  ANGULAR  MOTION 
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FIGURE  15.  ELEVATION  SYSTEM  BLOCK  DIAGRAM 


9.  Colculotion  of  Angie  £ 

The  geometry  of  angle  is  shown  in  Figure  16. 


GUN 


FIGURE  16.  GEOMETRY  OF  ANGLE  /3 


From  Figure  16, 


P 


sin 


O'  =  sin-1 


D  =  D' 

(75) 

+  e" 

(76) 

-i  t' 

i 

(77) 

1  D  3 

R  ■ 

(78) 

Now  from  trigonometry  of  oblique  triangles, 


3This  geometry  was  largely  developed  by  R.  E.  Yates. 


(79) 


/Ts  -  a)  (s  -  b)  (s  -  c) 

r  _  , - - 

V  s 

a  =  X 

b  =  R 

c  =  R' 

S  =  1/2  (a  +  b  +  c) 

Since  R  and  R'  are  fixed,  the  angle  /3  is  uniquely  determined  as  a  func¬ 
tion  of  X  . 

The  list  of  nominal  values  used  in  the  simulation  is  given  in  Table  I. 
TABLE  I.  LIST  OF  NOMINAL  VALUES  USED  IN  THE  SIMULATION 
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